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Abstract 

Powerful spectrum decision schemes enable cognitive radios (CRs) to find transmission opportunities in spectral 
resources allocated exclusively to the primary users. One of the key effecting factor on the CR network throughput is 
the spectrum sensing sequence used by each secondary user. In this paper, secondary users' throughput maximization 
through finding an appropriate sensing matrix (SM) is investigated. To this end, first the average throughput of 
the CR network is evaluated for a given SM. Then, an optimization problem based on the maximization of the 
network throughput is formulated in order to find the optimal SM. As the optimum solution is very complicated, 
to avoid its major challenges, three novel sub optimum solutions for finding an appropriate SM are proposed for 
various cases including perfect and non-perfect sensing. Despite of having less computational complexities as well 
as lower consumed energies, the proposed solutions perform quite well compared to the optimum solution (the 
optimum SM). The structure and performance of the proposed SM setting schemes are discussed in detail and a 
set of illustrative simulation results is presented to validate their efficiencies. 

Index Terms 

Cognitive radio, spectrum handover, maximum average throughput, sensing matrix. 
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OPPORTUNISTIC spectrum sharing has been developed through the new promising concept 
of Cognitive Radio (CR), in order to meet ever-growing spectrum demands for new wireless 
services. Conceptually, CR is an adaptive communication system which offers the promise of intelligent 
radios that can learn from and adapt to their environment [[Tl. The major issue in designing a cognitive 
radio network is to protect incumbent/primary users from potential interference problems while providing 
acceptable quality-of- service (QoS) levels for secondary users (i.e., unlicensed users). To this end, sensing 
capability is exploited in CRs which enable them to find some transmission opportunities called white 
spaces, i.e., temporarily-available spectrums which are not used by primary users (PUs). Limited number 
of possible observations and dynamic nature of observed signals lead to imperfect sensing which is usually 
described by false alarm and miss detection probabilities. The false alarm occurs when the PU is idle, 
but the Secondary User (SU) senses the channel as busy. While the miss detection is occurred when the 
SU senses an occupied channel as free. 

Average throughput of the SUs is one of the most important performance metrics which depends on 
the candidate primary channels for sensing and transmission, and it must be considered in designing 
appropriate sensing schemes. Generally, there exists more than one channel to be sensed by a CR. As a 
result, sensing schemes are commonly divided into two categories, i.e., wideband sensing and narrowband 
sensing. Sensing is wideband when multiple channels are sensed simultaneously. These multiple sensed 
channels can cover either the whole or a portion of the primary channels [|2l. On the other hand, when 
only one channel is sensed at a time, the sensing is narrowband. Ease of implementation, lower power 
consumption, and less computational complexity lead to great interest in narrowband sensing. When the 
narrowband sensing is used, an immediate question arises: which channel should be sensed first? In other 
words, to achieve the best possible performance, the channels have to be sensed in an appropriate order 
determined by sensing sequence (SS). 

In [|3]|, the problem of joint optimization of sensing and transmission is addressed. Specifically, Zhao et 
al. in [3J proposed a decentralized slotted CR MAC protocol to grasp the optimal policies for spectrum 
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sensing and access framework through a partially observable Markov decision process. Minimizing the 
overall system time of a SU, which contains the average waiting time and the extended data delivery time, 
through load balancing in probability-based and sensing-based spectrum decision schemes is investigated 
in dH. In [m, H, and [|71, the procedures to determine the optimal set of candidate channels for sensing 
are first discussed and then the maximization of the spectrum accessibility through optimal number of 
candidate channels are investigated. In [[8]|, |l9l, and [fTOll . the sequential channel sensing problems are 
formulated based on maximizing the throughput of the SUs. While in these works, the optimum sensing 
times have been studied, the effects of the sensing errors have not been addressed. Setting a SS by 
prioritizing the various channels can play a major role in finding a transmission opportunity or equivalently 
expected SU's throughput. Channel prioritization has been considered in [1 IJ in which an optimal channel 
sensing framework for a single-user case including the sensing order and the stopping rule has been 
proposed. In [11 J, it has been also assumed that the SUs are allowed to recall and guess. Recall means 
the ability to go back and access a previously sensed channel and guess means accessing a channel that 
has not been sensed yet. In ifTll and [[T3l . a stopping rule has been developed to determine when to stop 
the sequential sensing procedure and when to start secondary transmission. In [14], the optimal SS has 
been derived for channels with homogeneous capacities, and it was shown that the problem of finding 
the optimal sensing sequence for these channels is NP-hard. The authors in |[T5l have suggested a SS 
which sorts channels in descending order according to their idle probabilities. In [[T6l , finding the optimal 
SS sequence has been investigated for a single-user case with the aim to maximize the SU's throughput. 
The problem of finding optimal SS for two CR users has been addressed in y/7J, in which an exhaustive 
search has been applied in order to find the best sensing sequences for the users at the expense of a huge 
computational complexity. To reduce the complexity associated with the optimum solution, the authors of 
ifTTI have proposed two low-complexity suboptimal algorithms with the achieved throughput close to the 
maximum possible value. 

In this paper, the problem of selecting proper spectrum sensing sequences for a cognitive radio network 
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(CRN) with multiple users is addressed. Our objective is to maximize the average network throughput. 
First, we assume a perfect channel sensing (i.e., error-free sensing) and formulate an optimization problem 
on spectrum sensing sequences of the SUs based on maximizing the average throughput of the network. We 
discuss the conventional solution as well as its computational complexity. Due to massive computational 
burden of the conventional optimization algorithm, a novel algorithm, which finds near-optimal solution, is 
proposed. The proposed algorithm, called sensing matrix setting (SMS), provides short-term and mid-term 
fairness among the SUs and offers near-optimal solution with tolerable computational complexity as well 
as relatively low consumed energy. Then, we consider the impact of sensing error on the SMS algorithm, 
and propose modified version of SMS algorithm, called MSMS algorithm. In addition, for the multiple 
access among the SUs, we apply the conventional p-persistent MAC within the MSMS algorithm, and 
call the extended algorithm as PMSMS algorithm. Structure, performance, and related spectrum allocation 
processes for the proposed algorithms are discussed in detail. 

The rest of this paper is organized as follows. In Section |IIl we describe the CR network considered 
and the related assumptions. In Section [III the throughput of the CRN for a given sensing matrix is 
formulated, and the conventional approach to find the optimal SM as well as its computational complexity 
are discussed. In Section |IVl the structure, computational complexity, and consumed energy of the novel 
suboptimal SMS algorithm are described in detail. In Section |Vl the modified version of the SMS algorithm 
is introduced. The PMSMS algorithm is described in Section |VIl Numerical results are then presented in 
Section IVIIl which validate our analysis and verify the advantages of the proposed algorithms. Finally, 
the paper is concluded in Section IVIIIi 

II. System Model 

We assume a time slotted CRN with Ng secondary users which attempt to opportunistically transmit in 
the channels dedicated to the Np PUs. As in [[T3l . [fT6l . [fTTl . and [[TSl . the SUs are time synchronous in 
time-slots with other SUs and with the PUs. When a PU has no data for transmission, it does not use its 
time-slots; and hereby provides a transmission opportunity for the SUs. That is, at the beginning of each 
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time-slot, a channel can be established as occupied or vacant. In order to find the transmission opportunities 
appropriately and to protect the PUs from harmful interference, the sensing process is performed at the 
beginning of each time-slot. We assume that the SUs are equipped with simple transceivers, so they are 
able to sense only one channel at a time. The SUs always have packets to transmit, and as a consequence 
they will start transmission when an opportunity is found. Each SU senses the channels according to its 
SS sequentially, i.e., the SU senses the first channel from the top of its SS for a predetermined time 
duration (channel sensing time), and then senses the second channel if and only if the first channel is 
found busy. This procedure will continue until a transmission opportunity is found. Moreover, as [fT6ll , 
we assume that the SUs are not able to "recall" which means that they cannot re-sense and transmit on 
a previously sensed busy channel. 

The SU might stop its transmission in a channel and try to choose a new one due to the presence of 
the PU or the availability of a better channel with a more appropriate transmission condition. In order to 
switch to a new channel, which is called spectrum Hand-Over (HO), a secondary device needs a specific 
and constant time duration Tho to prepare its sensing circuitry for the next spectrum sensing. 

For the SU, each slot contains two phases: 1) sensing phase, and 2) transmission phase. The sensing 
phase contains several mini- slots of duration r (sensing time of each channel). Sensing is carried out 
by the SUs in the mini-slots, and once the transmission opportunity is found, the transmission phase 
will be started. This kind of access, i.e., listen-before-talk (LBT), is a common method in many wireless 
communication systems, for example see the quiet period in IEEE 802.22 standard [fT9ll . The sensing 
procedure is performed in an order based on the SS provided by the secondary network coordinator. 
The SUs do not have the adaptive modulation and coding (AMC) capability; so they transmit with a 
constant rate, R, during the transmission phase. We define a sensing matrix (SM) as a matrix with the 
dimensions of Ng x Np, in which the i-th row contains the SS dedicated to the i-th SU. Given the primary- 
free probabilities, i.e., Pq j , 1 < j < Np, and predetermined false alarm and detection probabilities, our 
objective is to find the optimal (or near-optimal) SS of each SU, i.e., the optimal SM, in order to maximize 
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the CRN throughput. 

Fig. \T\ demonstrates the slotted timing structure of a SU and its sensing process. For the example 
considered in Fig. [U after sensing {k — 1) occupied channels, the SU senses the A;-th channel free and 
then transmits data on that channel until the end of the slot. The wasted time length, i.e., the time allocated 
for sensing and HO in this process is equal to r + (A; — 1) (r + Tho)- Thus, the time left in the slot for 
transmission is T — r — (A; — 1) (r + Tho), where T is the time slot duration. Generally speaking, each slot 
is composed of a sensing phase with the maximum length of r + {Np — 1) (r + Tho) and a transmission 
phase with the minimum length of T — T — {Np — 1) {t + Tho) - Let us define the slot effectiveness, denoted 
by Cs, as the ratio of the transmission phase length to the slot length. Hence, if a secondary user starts 
transmitting on the A;-th channel of its SS, the slot effectiveness is: 

transmission time t + (k — 1) (t + Tho) 

= T- m — J— = 1 ^ — ^ 

time slot duration I 

and 

B, = fixe. = fl(l-I±(^^iMl±^) (2) 

where is the average throughput of the SU if the A;-th channel of the SS is sensed to be free and 
chosen for the transmission. From ([T]) by the increase of k, the channel effectiveness is reduced. 

in. Optimal Sensing Matrix 

In this section, we evaluate the CRN throughput for error-free sensing case and discuss about the 
optimal sensing sequences of the SUs (or equivalently SM) for the network throughput maximization. 
The optimal sensing sequence for the CRN containing just one SU is derived in [fT5l . For that case, all 
required was to sort the spectrums based on their primary-free probabilities (i.e., the absence probability 
of the PU). But in the CRN with multi users, the impact of collisions among the SUs' transmissions 
has to be taken into account. Assume that 5" denotes the sensing matrix which contains Ns rows and Np 
columns with the element Sij indicating that the i-th SU senses the Sij-th spectrum in its j-th mini-slots. 
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For network throughput evaluation, we note that for each spectrum Sij two possible cases might occur. 
First, the Sj j-th spectrum has been sensed by some SUs in their previous mini-slots (the mini-slots before 
j-th mini-slot). Li this case, regardless of the presence or the absence of the PU, the spectrum is occupied 
(as we assumed error-free sensing in this section). That is, the SU that senses this channel at the first time 
will transmit on the channel, as a result of perfect sensing, if the spectrum is idle. Second, the spectrum is 
sensed at the j-th mini-slot by i-th SU at the first time. So, the occupation probability of this channel only 
depends on the PU activity. If the Sij is sensed free, the i-th SU starts its transmission in this spectrum 
for the rest of the time slot with a constant rate R. From the above discussion, repetition of a spectrum in 
the sensing matrix when assuming error-free sensing does not offer any benefit to the CRN throughput. 

Fig. [2] shows the structure of the SM. In this Fig., demonstrates the spectrums allocated prior to 
the j-th mini-slot. The array Y contains the spectrums dedicated at the j-th mini-slot to the SUs prior to 
the user i. Let Ag-J^ indicate the presence or absence of the spectrum Sij in Zj^i, 



By the above definition, when not taking into account the impact of collision caused by multi SU 
transmissions at the same spectrum due to simultaneously finding the spectrum free at the j-th mini-slot, 
the spectrum Sij can be efficiently used by the i-th SU with the probability of Pq ^. ^^As- J^ , where Po,Si j is 
the absence probability of the Sij-th PU. To consider the impact of the collision on the network throughput 
as discussed above, we define the operator © as follows. 




(3) 
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The operator © is used to model the possible collision due to multi SUs finding the same spectrum 
free at the j-th mini-slot. As stated before, each channel in each time-slot has a contribution in the whole 
throughput if and only if it is sensed only by one SU (i.e., assigned to one SU sensing sequence) because 
of error-free sensing assumption. 

By the above definition of the operator ©, the average throughput of the CRN is easily computed as 
follows, 

{Po,s,..K\ © ^o,s...<, © ■ ■ ■ © Po,s^^.,A^l^,) i?2 + ■ ■ ■ + (6) 

where Bj is defined in Q and Af."^ = 1 ,Vz ^ A'p. Now, the optimal SM is found by solving the following 
optimization problem, 

S* = argmax Q (7) 

Sl.l,Sl,2,...,SNs, Np 

Based on the derived optimization problem, we can find the optimal SM by exploiting exhaustive 
search. Assume that the computational complexity of computing © for a given SM, S, is in 0{l). 
Then, the computational complexity of finding the optimal SM is in O (^Np^^^^^). Since the expression 
describing the performance metric (the CRN throughput) is complicated in general, there is no much 
room for solving dV]) through classical optimization procedures. On the other hand, solving (|7]) through 
the exhaustive search makes no guarantee for fairness among the SUs. In addition, it results in massive 
computational burden which is not scalable regarding to both Np and Ns. All these facts make a strong 
motivation and interest in developing an appropriate suboptimal solution for the problem formulated in 
d?]). In the following sections, we propose suboptimum solutions for the SM for three different cases. 
Advantages of proposed algorithms are threefold: First, it offers low computational complexity. Second, it 
provides fairness among the SUs. Finally, its consumed sensing energy to find a transmission opportunity 
is much less than the exhaustive search. 
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IV. SMS Algorithm 

A. Structure of the SMS Algorithm 

The proposed algorithm, designed for error- free sensing case, is composed of Np rounds. In the k- 
th round, the coordinator determines the k-th column of the SM, i.e., Si^k, 1 < ^ < Ng. As mentioned 
before, repeating a spectrum in the SM for more than one times either in the same mini-slot or in different 
mini-slots does not have any benefits on the network throughput. During each round and for each SU, 
the coordinator assigns a reward to each candidate channel^ to be possibly allocated to the SS of the SU 
at that round and then adopts the channel with the maximum reward. That is, at the round m, for each 
secondary user k and for each unassigned channel i, we define Cf'^ (m) as the reward of the channel 
i if selected as the m-th component of the k-th SU's sensing sequence. This reward is set equal to the 
contribution of the k-th SU to the network throughput if the i-th channel is selected, as will be described 
latter. Then, the channel with the maximum reward is selected. 

We denote the set of all assigned channels to the sensing matrix by A. At the beginning, we have 
y4 = and 5 = , where S is the sensing matrix and denotes empty matrix. We also denote the set of 
all channels by N, where N = {1, 2, . . . , Np}. The process is as follows: 

Round- 1 

For this round, first the coordinator assigns a spectrum to the SS of the first SU at its first mini-slot. The 
coordinator must adopt si i from the unassigned channels, i.e., A = N\A = N. We have (1) = Po^iBi, 
where Bi is defined in Q. The coordinator selects a channel with the highest reward for the Si i. That 
is, the first channel to be sensed by the first SU is, 

si,i = argmax Cf-* (1) (8) 
ieA 

After si 1 is determined, A and A are respectively updated to A = {si i} and A = N\A. This procedure 

'the channel that has not been assigned to the SS of any user previously 
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is repeated for each SU; so for the i-th user in the first round, we have, 

S£^i = argmax cf'' (1) (9) 

where cf^ (1) = Po,iBi. 
Round-m 

At the m-th round, for each SU, the coordinator similarly assigns a reward to each left spectrums and 
allocates the best spectrum, which has the maximum reward, to the SU. If the coordinator chooses the 
j-th channel for the m-th sensing mini-slot of the i-th SU, the following reward will be gained by the 
user. 

Gf (m) = (1 - ^o,s,J j Po,B^ (10) 

Hence, the coordinator determines the m-th element of the SS of the £-th SU as, 

S£,TO = arg max Cf ^ (m) (11) 

At this round, first it must be determined from which SU the procedure should be started. In order to 
achieve an acceptable level of fairness among the SUs, the algorithm starts with a SU that has gained 
the lowest cumulative rewards during the previous (m — 1) rounds (previous (m — 1) mini-slots). The 
cumulative reward gained during the previous (m — 1) mini-slots is calculated for the i-th SU as, 

m—l 

E (^) = ^£ (1) + (2) + ■ ■ ■ - 1) (12) 

k=l 

where G^se,^ (k) is defined in (fTOl) . 

This process continues until \A\ = Np or equivalently A = 0. At the end of the process, the elements of 
S without any assigned spectrum are replaced by zero, which indicates that the sensing is not performed 
for those elements. Since each channel is sensed only once in the proposed algorithm, the energy consumed 
by the SMS algorithm equals to NpEc at the worst case, and it does not increase by A^^. 

In order to have mid-term fairness, at the beginning of second time-slot, the process starts with the 
second SU and the first element of the sensing sequence of this user is determined and then the procedure 
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is continued by selecting the first element of the third user, and at the last the first element of the first 
user is selected. The other elements are determined as described above. This cyclic ordering is continued 
in the following time-slots, i.e., at the beginning of m-th run of the SMS procedure (m-th time-slot) the 
process starts with selecting the first element of the sensing sequence of the k-th SU, where 

k = mod {m, Np) (13) 

These procedures are summarized in Algorithm [U 

B. Computational Complexity 

As we stated before, the computational complexity of finding the optimal SM is in order of O [Np^''^^"), 
while it is in order of O (1) for our proposed method. In the SMS algorithm, a channel will be assigned 
to the SM if it offers the highest reward, defined in ([10), among the left channels. From (flOl) . it can be 
easily shown that G/''^ (m) > Gk^''^ (m) if Pqj > Po,fc □• So for the the error-free case, the information 
required to determine the SM is the primary-free probabilities of the channels. 

C. Averaged Consumed Energy for Finding a Transmission Opportunity 

Let Ec (r) and Ec (rho) denote the consumed energy for sensing of each primary channel and the con- 
sumed energy for each HO, respectively. Hence, the average consumed energy for finding a transmission 
opportunity can be calculated as, 

{Ns + gi+g2 + --- + gNj E, (r) + {gi + g2 + ■ ■ ■ + ^ivj E, (rho) (14) 

where gi denotes the average number of HOs required by the i-th SU to find an idle channel. 

The processes of channel sensing and signal transmission consume more energy compared to the HO. 
Therefore, it is rational to ignore the second term {gi + g2 + ■ ■ ■ + gNs) Ec (rho) in (fT4l) compared to the 
first one. 

^The reason for defining the award as in dlOl l is that it can be easily modified to the non-error-free sensing case and also for the case of 
considering different MAC schemes. 
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Algorithm 1 sensing matrix setting algorithm for perfect sensing case(the SMS algorithm) 
Initialization: 5 = 0, N = {1, 2, . . . , Np}, A = ^,A = N\A, and m = 1 

H i— sort the SUs based on their numbers, sequentially. 

for m = 1 to (the maximum number of the SUs' mini-slots) do 

for i = (the first element of H) to (the last element of H) do 

for z = 1 to A'p do 

if A ^ then 

Compute Cf^ (m) as in (fTOl) 

Assign Si^rn ^ argmax Cf^ (m) 

else 

se,m ^ 

end if 

end for 

if the spectrum i is assigned to m then 

Add {i} to A 

Update A 
end if 
end for 

H ^ sort the SUs based on their cumulative rewards computed using (fT2)) . 
end for 

Return S as the suboptimal sensing matrix. 
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To evaluate the average number of HOs of the i-th SU, (ji, we consider two following cases: 1) The SU 
searches among the channels, finds a transmission opportunity, and then transmits, 2) The SU searches 
among the available channels, but does not find any free channel. Then, gi can be easily calculated as, 

Np-l Np 

= Pm,,Po,s,. + 2Pi,,,,,Pi,,,,,Po,s,3 + ■ ■ ■ + (iV, - 1) Po,.,,^ n Pi,.,, + iV, n Pi,s.„ (15) 

3=1 j=l 
. ' 

Tl 

where the term Tl represents the average sensed channels by the i-th SU until the user finds a transmission 
opportunity, and the last term demonstrates the case that the SU senses all channels busy, and therefore 
it does not transmit on any channels assigned to its SS. By substituting (fT5l) in (fT4l) . the total average 
consumed energy for the exhaustive search method is derived as follows: 

Ns Ns / Np / k-1 \ \ Ns Np 

X; 9^Ec = £ 1 + X] (/^ - 1) n ^M.. ^0.^.. ]+EcJ^N,l[ Pi,,,, (16) 

i=l i=l \ k=2 \ j=l J J 1=1 j=l 

On the other hand, for the SMS algorithm, analytically deriving the average consumed energy is 
complicated. Hence, we only focus on two special extreme case, i.e., the maximum and minimum 
consumed energies. For the worst case, which consumes the maximum energy, all the channels appeared 
in the SM are sensed. In this case, the consumed energy equals to NpEdr). For the best case, which 
consumes the minimum energy, each SU finds the first channel of its SS free and does not need to sense 
the rest. In this case, the consumed energy is equal to min (A^p, A^.) Ec (r). It is worth noting that for 
Pqj = 1 , for all j, the sequential sensing scheme forces the SUs to continue searching among all channels 
in their sensing sequences, which is equivalent to the worst case, and similarly Pqj = , for all j is 
equivalent to the best case with minimum consumed energy. If we compute the average consumed energy 
of the optimum solution given in (fT6l) for theses two cases, the maximum and minimum consumed 
energies will be equal to EcNg (1 + Np) and EqNs, respectively, which are higher than those of the SMS 
algorithm. 

For more elaboration, we study an especial case where all channels have the same primary-free 
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probabilities, i.e., Pqj = P , Pij = 1 — P , for all j. Then, we can simplify (fT6l ) as follows, 

Ns / Ns / Nj, \ 

Ec {r)Y^h = Ec{t)[NsP + pY,[Y. ((^ " 1) " P)'') +NsN,{l - Pf^ 



1=1 \ 1=1 \k=2 



1 - P + (-1 + P - PNp) (1 - P)^" 



Ec (r) Ns\P+ ^ ^ ' ^ ^ ' ^ ^-'''^ ^ ' + N,{\ - Pf- 1 (17) 



It is worth noting that (flTI) is a decreasing function of P. Hence, the minimum and maximum values 
of consumed energy are related to the cases P = and P = 1, as discussed above. Moreover, as will be 
shown in the numerical result section (Fig. [5]), the consumed energy associated with the optimal SM is 
higher than the consumed energy for the matrix obtained by the SMS algorithm for all values of P. 

In the following, the impacts of the sensing errors are investigated. In general, the sensing error manifests 
itself in two forms: false alarm and miss-detection. In the SMS algorithm proposed, a channel is allocated 
to only one SS, and thus the SUs have no common channels in their sequences. Although this approach 
performs well when there is no sensing error, but in the case of non-perfect sensing, this method is not 
efficient; since by a false alarm made by a SU in a sensed channel, a transmission opportunity is lost 
for this channel by all SUs. Therefore, it seems that the coordinator has to repeat spectrums in the S in 
order to increase the possibility of exploiting all opportunities and thus to increase the spectral efficiency. 
On the other hand, allocating a channel to the sensing sequences of multiple SUs increases the average 
number of sensed channels and thus raises the average sensing energy consumption. Moreover, due to 
miss-detection, it is possible for a SU to mistakenly transmit on a channel which is already used by 
another SU or PU, and therefore some collisions might occur. Hence, there is a trade-off between the 
average achievable throughput, energy consumption, and the level of collision in the CRN which must be 
addressed in an extension of the SMS algorithm. 

To modify the SMS algorithm, we must consider the impact of sensing error probabilities on the reward 
function. Moreover, the channels occupation probabilities will be different at the beginning of the various 
mini-slots, which must be reflected in the reward function. Finally, because of repetition of each channel 
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in the SM, the stopping rule, which was \A\ = Np, must be modified. 

V. MSMS Algorithm 

Since in the SMS algorithm, the sensing sequences of the SUs have no common channels, the occupan- 
cies of channels at the beginning of each mini-slot only depend on the PUs' activities. But if the channels 
are allowed to be repeated in multiple rows or columns of the SM, the occupancy of a channel can be due 
to the presence of either the PU or a SU. To extend the SMS algorithm, first the occupation probability 
of the j-th channel, i.e., qij, has to be determined. 

It is worth noting that, since all the SUs use the same sensing schemes with the same sensing time 
lengths, they all have the same probabilities of false alarm and miss -detection. Thus, we have, 

Pfa,l = Pfa,2 = ■ ■■ = Pfa,Np = Pfa 

(18) 

Pd,l = Pd,2 = ■ ■ ■ = Pd,Np = Pd 

In order to reflect the impact of sensing error on the proposed algorithm, three possible cases must be 
considered when the coordinator tends to adopt the /-th channel as Sij: 

• The /-th channel has not yet been allocated to any sensing sequences, i.e., / ^ (matrix Z is 
defined in section [III] and Fig. |2l). 

• The /-th channel has been adopted at least once for sensing at the previous mini-slots, i.e., / G Zj_i. 

• The /-th channel has been allocated simultaneously to multiple users at the j-th mini-slot (vector Y 
shown in Fig. 2). 

Fig. [2] depicts these cases graphically. Suppose that j-th component of the SS of the i-th SU, i.e., Sij, 
is to be selected by the coordinator. So the reward gained by adopting the fc-th channel as Sij is to be 
determined. Considering the definition of the matrix Z and the vector Y in Section |lll] and also in Fig. |2l 
if n elements of Z are equal to k, this will indicate that the fc-th channel has been sensed at most by n 
SUs during previous mini-slots. Also, if two or more elements of Y are equal to k, then the k-th channel 
will be sensed by two or more SUs during the j-th mini-slot. When channels are allowed to be sensed 
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by multiple SUs simultaneously, an appropriate MAC protocol can be used to regulate the access of the 
SUs to transmission opportunities. As a first step, we assume that a SU starts transmitting when it finds a 
transmission opportunity. Applying more appropriate MAC protocols to decrease the collision probability 
among the SUs will be considered later. For the mentioned transmission policy, if the k-th channel belongs 
to Y and it is also adopted as Sij by the coordinator, a collision may occur and the reward may be zero. 
Without loss of generality, it is assumed that the coordinator starts the allocation process for the j-th 
mini- slot from the top of the j-th column of matrix S. 

Given Zj^i, we denote the occupation probability of the k-th channel at the beginning of the j-th 
mini- slot as q[^l. Then, we easily obtain, 

q[',l = Pi,k + 9^,,,z,., (19) 

where 6i^k,Zj.i represents the probability of several sensing of and possibly transmitting on the k-th 
channel in the first (j — 1) mini-slots and is easily computed as, 

9iAz,., = E St^'^-' (1 - Po,kPfan (20) 

n=l 

where \Zj^i\ = Ns (j — 1) is the number of elements of Zj^i and 6n^^'^ is defined as, 

{1 : if spectrum # A; is in for n times 
(21) 
: if spectrum is not in Zj_i for n times 

As in the SMS algorithm, the process starts with 5 = and at the first step si^i is selected for the SS 
of the first SU by the coordinator. As before, for channel i E N, (1) denotes the reward contributed 
by the A;-th SU to the overall throughput of the secondary network when the i-th channel is allocated as 
the first element of its sensing sequence. 

Round- 1 

At the first round of the MSMS algorithm, cf^ (1) = (1 - Pfa) Bi in which Bi is defined in ©. 
ql^^ = 1 — q^j, where q^J is given in (fT9l ). Therefore, si i is determined as, 

Si^i = argmax G.^-* (1) (22) 

iGN 
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If the i-th spectrum is adopted as the reward added to the system can be calculated as 

\y\ 

gT (1) = E ^'^<&}Pf^'' (1 - Pf^) Bi (23) 

n=0 

in which 

{1 : if spectrum #i is in Y for n times 
(24) 
: if spectrum #i is not in Y for n times 

Finally, the first channel of the SS of the Z-th SU is selected according to: 

s,,i = argmax Cf (1) (25) 

Round-m 

The l-th SU gains a reward by adopting the j-th channel as the m-th element of its sensing sequence 
provided that the user has not detected a transmission opportunity in its previous sensed channels. Note 
that besides finding a truly free channel, the SU may mistakenly sense an occupied channel as free due 
to miss detection and does not continue the sensing procedure. Therefore, the reward gained by the l-th 
SU is: 



Cf (m) 



m— 1 



K 5: '/ 



\Y\ 
n=0 



C2 



where CI indicates the probability of requiring (m — 1) HO, and C2 represents the throughput contribution 
of j-th channel if selected at the m-th mini-slots of the l-th SU for the transmission. 
Thus the m-th element of the l-th sensing sequence is similarly determined as, 

Si^rn = arg max cf^ (m) (27) 

Similar to the SMS scheme, in the MSMS algorithm, at the round-m (m G {2, 3, . . . , Np}) the coordi- 
nator starts with the SU that has gained less cumulative rewards in its (m — 1) previous mini-slots. The 
cumulative rewards of the j-th SU at its (m — 1) previous mini-slots can be computed as, 

Gif , (1) + Gif^ (2) + • • • G^l_^ (m - 1) (28) 
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where G^^^ (z) for 1 ^ z ^ m — 1 is calculated as (|26l ). Hereby, a certain level of fairness is ensured 
among the SUs. 

The stopping rule of the MSMS algorithm is different from that of the SMS algorithm. For the MSMS 
algorithm, two possible rules can be exploited. First, there exist no constraint on the number of times that 
each channel can be used as the elements of the SM. For this case, the process is stopped when all the 
elements of the SM have been selected. Second, the number of times that each channel is appeared in 
the SM is limited. While the first rule leads to the maximum average throughput which can be achieved 
by the MSMS algorithm, the second rule is more rational and practical. The probability of a channel 
erroneously sensed as busy exponentially decreases by the number of times that the channel is sensed. 
As a result, we use the second stopping rule. In the numerical result part, we limit the number of times 
that each channel is appeared in the SM to 3. That is, the coordinator assigns each channel, if necessary, 
at most three times in the SM. 

In order to have further mid-term fairness among the SUs, the same idea as applied to the SMS algorithm 
is exploited, i.e., at the beginning of m-th run of the MSMS procedure (m-th time slot), the process starts 
with the A;-th SU as specified in (fT3l) . The procedures of MSMS algorithm are summarized in Algorithm 
121 

VI. PMSMS Algorithm 

Regardless of how the SM is created, it is possible for a channel to be assigned to the several SUs in 
the same mini-slot. In this case, various conventional MAC algorithms can be exploited to increase the 
transmission chance on this channel. In this section, we utilize the well-known p-persistent MAC (PMAC) 
protocol in the MSMS algorithm and develop PMSMS algorithm. In this algorithm, in each mini-slot the 
SUs sense the assigned channels with the probability of p. In order to have a synchronous sensing scheme 
for all SUs, the SU will be idle for r seconds(mini-slot time duration) if its MAC protocol does not allow 
it to sense the channel. 

The stopping rule as well as fairness establishment techniques are similar to the MSMS algorithm. 
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Algorithm 2 sensing matrix setting algorithm for non-perfect sensing case(the MSMS algorithm) 
Initialization: 5 = 0, N = {1, 2, . . . , Np}, A = N,m = l, and RN = {0,0,. .. , 0}^^^^ 

H ^ sort the SUs based on their numbers, sequentially. 

for m = 1 to (the maximum number of the SUs' mini-slots) do 

for i = (the first element of H) to (the last element of H) do 

for i = 1 to Np do 

if A ^ then 

Compute q[f using and ^ 

Compute Cf ■* (m) as in (|26|) 

Assign Si^m ^ arg max Cf ^ (m) 

else 

se,m ^ 
end if 
end for 

if the spectrum i is assigned to si^m then 

RN [i] ^ {RN + 
end if 

if RN \i] = 3 then 

A = n\ {i} 
end if 
end for 

H ^ sort the SUs based on their cumulative rewards computed using (|28]) . 
end for 

Return 5" as the suboptimal sensing matrix. 
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Considering PMAC, there are two cases that a free channel is not used by a SU. First case is due to 
the false alarm, and the second is due to the presence of PMAC protocol. In the latter case, the channel 
is not sensed with the probability of (1 — p). Considering these two cases easily leads to the following 



modification of 9i^k,z^i (defined in (|20l)). 



^l,fc,Zj_i 



n=l 

\Zj-i\ 



1-E 



n 



t=o 



\ 



\ 



n—t 



^ n=l 



, <p < 1 



p = 1 



(29) 



Then, the channel occupation probability is obtained by substituting (|29] ) in (fT9l ). The generalized reward 
of assigning the j-th primary channel to the SS of the i-th SU at the m-th mini-slot is simply calculated 



as. 



' m—l 



\Y\ 



Gf (m) 



(30) 



i=l 



n=0 



where 



il-Pfa)pE 
k=0 



n 



Pfa'il-Pfar-'il-vT-^ 0<p<l 



(31) 



(1 - Pfa) Pfa^ , P=l 

and finally the coordinator adopts a channel with a highest reward for the m-th mini-slot of l-th SU as 
follows, 

se,m = argmax Cf-* (m) (32) 

The SUs sense the assigned channels with more probability as p increases, which can increase the 
chance of finding a transmission opportunity, at the expense of raising the level of contention among the 
SUs. Therefore, there is a tradeoff on the value of p, which will be discussed in the next section. 



VII. Numerical Results 

In this section, the performance of the proposed allocation schemes is evaluated by simulation consid- 
ering the effect of different parameters. Moreover, advantages of exploiting the proposed algorithms are 
demonstrated through exhaustive simulations. 
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The simulation parameters are given in Table HI The values of SNR and sampling frequency which 
are used by the energy detector are adopted from IfTSl . The value of sensing time of each channel, r, is 
selected such that the false alarm and detection probabilities meet the constraints imposed by the IEEE 
802.22 standard |fT9l . Each SU senses the channels according to its sensing sequence, each for r seconds, 
until a free channel is found. Then, the SU transmits on this channel for the rest of the time slot. The 
average normalized CRN throughput has been evaluated by simulating the scenario for 100 time slots. 

Fig. [3] validates our analysis and depicts the average throughput (normalized to R) of the SUs versus the 
sensing time, for the optimal SM and the SM obtained based on SMS algorithm, for a error-free sensing 
case. For the optimal SM, both the theoretical and simulation results are provided. As it can be realized, 
the throughput linearly decreases by the sensing time. This is due to the fact that for the error-free case, 
there exists no error in the detection scheme and thus while the increase of sensing time does not have any 
positive impact on the correct detection, it linearly reduces the transmission time, as can be inferred from 
Fig. [3] also verifies near-optimality of our proposed algorithm; while it imposes much less complexity 
burden than the optimum scheme. The relative difference between the average throughput obtained by the 
SMS algorithm and that obtained by the exhaustive search method is negligible and about 0.81%. 

Fig. m compares the average throughputs of different SUs for a CRN with three SUs, again for a error- 
free sensing scheme. In this example, it is assumed that the number of primary channels, Np, is 5. The 
maximum relative difference between the SUs throughputs is 1.84% which confirms the fairness among 
the SUs when using the proposed scheme. It is expected by running the simulation for more than 100 
times, the difference among the SUs' throughputs disappears. 

Fig. [5] demonstrates the energy efficiency of our proposed algorithm. This Fig. shows the average 
consumed energy versus the primary user occupation probability. The SUs consumes less average energy 
to find a transmission opportunity when sensing the channel based on the SM obtained by our proposed 
method compared to the optimal SM. It is worth noting that in both schemes, the consumed energies of 
the SUs increase when the PUs' absence probability decreases; as the SUs have to sense more channels 
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to find free ones. 

Fig. |6] compares the average throughputs of various spectrum allocation schemes proposed in this paper. 
For the two MSMS and PMSMS schemes, a practical scenario with activity detection errors has been 
assumed. In general, as the sensing time increases, the detector senses the channels more accurately and 
finds more transmission opportunities. However, by the increase of the sensing time, less time remains 
for the transmission. Hence, there exist a tradeoff between average throughput and detector accuracy. As 
it is seen in Fig. [6l first the SUs' throughput increases by r (due to an accurate sensing); then after 
an optimum point, where Pmd and Pja are in acceptable levels, the throughput starts decreasing due to 
the reduction of the time left for the transmission. For a sensing time greater than a specific amount 
(optimum value), the false alarm and miss-detection probabilities of the detector becomes negligible, and 
the allocation procedure of the MSMS algorithm as well as its performance will be similar to those of the 
SMS algorithm, for which a error-free sensing has been assumed. In the MSMS scheme, the SUs sense 
channels with the probability of (1 — Pfa), and thus some transmission opportunities are lost as a result 
of false alarm. This is the reason that the average throughput of the SUs obtained by the MSMS algorithm 
is less than that of the SMS algorithm in which Pfa is assumed to be zero. In the PMSMS algorithm, 
the applied p-persistent MAC protocol leads to loss transmission opportunities, and thus to less average 
throughput compared to the MSMS when the number of SUs is not too high. 

Fig. |7] demonstrates the advantages of the exploited PMAC for the case that the number of primary 
channels are less than the number of SUs. This Fig. shows the average SUs' throughput versus the 
probability of sensing an assigned channel (i.e., p in MAC protocol). Note that the performance for p = 1 
is the same as that of the MSMS. As can be realized, for A^^ = 8 and Np = 5, the exploited PMAC 
protocol can increase the chance of transmission on the channels by reducing the contention level among 
the SUs. Therefore, the PMSMS scheme can offer higher throughput for the CRN than the MSMS scheme 
provided that the coordinator appropriately selects the value of p, which for the example considered it 
must be larger than 0.165. Interestingly, for p = 0.46, the improvement in the throughput when using the 
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PMSMS scheme is about 48.8% compared to the MSMS scheme. 

VIII. Conclusion 

In this paper, the average throughput of a cognitive radio network (CRN) for a given sensing matrix 
(SM) has been derived, and an optimization problem has been formulated to find the optimal SM. In 
order to mitigate the challenges associated with the optimal solution, three novel centralized suboptimal 
algorithms have been proposed. More specifically, the SMS and MSMS schemes are proposed for error- 
free and non-perfect sensing cases, respectively, and then the PMSMS algorithm is developed by applying 
the conventional p-persistent MAC protocol in the MSMS scheme to strengthen the multiple access 
capability of the CRN. Besides offering throughput close to the maximum achievable one, the benefits 
of these proposed schemes are threefold. In addition to comparatively low computational complexities, 
they provide an acceptable level of fairness among secondary users. Further, they offer lower consumed 
energies compared to the optimum solution. The performance of the proposed schemes has been evaluated, 
and their efficiencies have been demonstrated through theoretical analysis as well as exhaustive simulation 
results. 
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TABLE I 

Simulation Parameters 



Parameter 



Description 



omax 



Is 

T 

Tho 
Np 



Value 



Minimum allowable detection probability 0.9 

Maximum allowable false alarm probability 0.1 

Receiver sampling frequency 6 MHz 

Time-slot duration 200 ms 

Required time for handover 0.1 ms 

Number of primary users 5 

Number of primary users 3 



time-slot duration, T 



Sensing of k-th 
channel 



channel sensing 
data transmission 



Fig. 1. General timing structure of our system model. 




Fig. 3. Average throughput versus sensing time for various sensing matrix selection schemes 
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Fig. 4. Average throughput versus sensing time for different secondary users 
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Fig. 5. Average consumed energy versus the primary user absence probability 
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Fig. 6. Average throughput of SUs versus sensing time for various proposed schemes. 
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Fig. 7. Advantages of exploiting PMAC. 



